Zero modes of tight binding electrons on the honeycomb lattice 
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Tight binding electrons on the honeycomb lattice are studied where nearest neighbor hoppings in 
the three directions are ta,tb and tc, respectively. For the isotropic case, namely for ta = tt — tc, 
two zero modes exist where the energy dispersions at the vanishing points are linear in momentum 
k. Positions of zero modes move in the momentum space as ta,tb and tc are varied. It is shown 
that zero modes exist if|^— 1 < |^ < 7*-+l- The density of states near a zero mode is 
proportional to \E\ but it is propotional to y'|i5 at the boundary of this condition 



PACS numbers: 81.05.Uw, 71.20.-b, 73.22.-f, 73.43.Cd 
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The integer quantum Hall effect has been observed in 
graphene^'^ when the carriers are changed by the gate 
voltage. The quantization of the Hall effect is observed 

2 

as (Jxy — 2n^ with n = ±1,±3---, where the fac- 
tor 2 comes from the spin degrees of freedom. These 
quantum numbers are unusual, since in a usual case 
n — 0,±1,±2,- • •. This unusual quantum Hall effect was 
discussed in terms of relativistic Dirac theorji^. However 
it is more natural to be explained by the realization of 
the quantum Hall effect in periodic systems^ in the pres- 
ence of zero modest. We will call zero modes instead 
of massless Dirac excitations in this paper because we do 
not consider relativistic particles. The energy spectrum 
and the density of states of the honeycomb lattice near 
half filling and in zero or small magnetic field are similar 
to these in the square lattice near half filling in a very 
strong magnetic field about half flux quantum per each 
unit cell^. 

At zero carrier concentration (i.e. half-filled elec- 
trons), the resistivity pxx is close to the quantum value 
/i/(4e^) — 6.45A:il independent of temperaturei, which 
has been also attributed to the zero modeaiSii. 

The existence of zero modes has also been proposed for 
the quasi- two-dimensional organic conductor a-(BEDT- 
TTF)2l3. The conductivity under pressure is almost con- 
stant in a wide range of temperatur^. Pertinent numer- 
ical computations performed by Kobayashi et aliS found 
that, for certain range of parameters, the Fermi surfaces 
become points and the density of states is proportional 
to energy at 3/4 filling of electrons. The existence of 
zero modes was also confirmed by the band structure 
calculationifliii. The unit cell for the model of a-(BEDT- 
TTF)2l3 has four non-equivalent sites. Katayama et al.^^ 
studied simpler model with two sites in the unit cell and 
they obtained a condition for zero modes. 

In this Letter we study a tight binding model on the 
honeycomb lattice and obtain the condition of ta, tb and 
tc for the existence of zero modes. 

Unit cell of the honeycomb lattice contains two sub- 
lattices as shown in Fig. ^. The Bravais lattice is a 
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FIG. 1: (color online) a. honeycomb lattice. Unit vectors are 
vi = (^, -^) and V2 = (^, ^). Three nearest neighbor 
hoppings are ta, tb and tc b. The red hexagon is a Brillouin 
zone for the honeycomb lattice. The reciprocal lattice vectors 
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where a is a distance between nearest sites. We consider 
only nearest neighbor hoppings. There are three nearest 
neighbors for each site, ta, tb and tc as shown in Fig^ We 
study the generalized honeycomb lattice model where ta, 
th and tc are not necessarily equal. Under uniaxial pres- 
sure, ta, tf) and tc have different values for each other. For 
example, ta> t^ ~ tc is expected, if the uniaxial pressure 
along the x direction is applied. The Hamiltonian for the 
generalized honeycomb lattice is given by 



n^Y^ "y- ta{alj^„^ + h.c.) 
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Using the Fourier transform 
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where x = (a,0), we obtain 
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The energy is given by 



^k — ^n + th + ^r 
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If we perform a translation in the momentum space 



(Kxi f^y) ^ (f^x I 'q"'' "'jy/' 
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and a replacement ia ^ — ^a simultaneously, we get the 
same Ck- Therefore we can take ta > without loss of 
generality. In a similar way one can take ib > and 
tc > without loss of generality by taking a translation 
in the momentum space, 
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[kx.ky) -^ (kx + -TT, ky ± -T-T^) 



The reciprocal lattice vectors are 
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as shown in Fig. ^p. Let us write 

k = fciGi + fc2G2, 



where 
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The energy is 

<^l = tl+tb+tl + '^tJb cos(27rfci) 

+ 2tatc cos(27rA:2) + 2tbtc cos(27r(-fci + ^2))- (15) 
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FIG. 2: Graphical explanation for appearance of zero modes, 
(a) If ta, tb and tc do not form a triangle, there are no zero 
modes and gaps at _E = are open, (b) Zero modes exist 
when ta, tt and tc form a triangle. Angles fci and fc2 are 
shown. 




FIG. 3: Zero modes exist in the filled region. 



The minimum of |ek| is obtained as follows. Consider 
the quadrangle ABCD in Fig. |2t- Then we have 



\DA\^ = \AB + BC + CD\ 



tl + tl + tl 



2tatb COS /3 — 2thtc COS 7 — 2tatc cos(7r — /3 — 7) > 0. 

(16) 



Put 



then 



/3 = TT - 27rfci 

7 = 7r + 27r(fci - k2), 



(17) 
(18) 



tl+tl+tl + 2tatbCOs{27Tki) 

-2tatc cos(27rfc2) + 2tbtc cos(27r(A:i - fe)) > (19) 
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FIG. 4: (color online) Energy dispersion for the isotropic case 
{ta = tb = tc = 1) ■ (a) 3D plot and (b) contour plot. 



The equality is satisfied when ta, tb and tc form a triangle 
which can be seen in Fig. [2l3., i.e., 



cos(27rfci) 
cos(27rfc2) 
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2tata 
+ 2 + 2 ^2 



cos(27r(fci - k2)) - 2t,t, 
The triangle can be formed when 



(20) 
(21) 
(22) 



N 

\ta\ 



-1 



< J_^ < 
- \ta\ - 



N 
\ta\ 



1 



(23) 



is satisfied. See Fig. |31 

In the isotropic case where ta — tf, ^ tc, zero modes are 
at {ki,k2) = ±(|,|), (fci,fc2) = ±(f.i)> and (A;i,fc2) = 
±(—■5-, 4) i-S. the corners of the first Brillouin zone, 



3' 3 



_ _i_(2TT 2V3tt\ 



(kx,ky) — ±( — , (j ), [kx,ky) — ±( — , g ) and 

ik^,ky) = ±(0, ^^). See Fig.H The density of states 
is plotted in Fig. |5^. 



FIG. 5: Density of states of the electrons on a generalized 
honeycomb lattice. 



If the parameters are in the boundary as seen in Fig. 
two zero modes merge into a confluent point. For ex- 
ample, Ef; = at confluent point (^11^2) = (0' 1/2) for 
ta — tb ^ I, tc — 2 (Fig.lHJ. Near this point ek is written 



ek oc ±Jciiki - kl)* + C2(fc2 - k^y 



(24) 



where ci and C2 are constants. In this case the density 
of states near E = becomes 



N{E) ex 



(25) 



See Fig. |Slb), while N{E) oc |-E| in the case of two zero 
modes (Fig. Et)- When the inequality Ea. H23|l is not 
satisfied, a finite gap opens at £' = as shown in Fig. [3;. 

In conclusion, we have studied the energy of tight bind- 
ing electrons in the generalized honeycomb lattice and 
found the condition for the existence of zero modes. The 
zero modes exist at the corners of the hexagonal first 
Brillouin zone for the usual honeycomb lattice. Two zero 
modes moved to become a confiuent point at the critical 
values of parameters ta, tb and tc, where ta, tb and tc 
stop to form a triangle. 
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FIG. 6: (color online) 3D plot (a) and the contour plot (b) 
of the energy of the generalized honeycomb lattice with ta = 
tb = I, tc = 2. 
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